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I .  INTRODUCTION 


This  study  vas  performed  to  iavestigate  and  document  the  effectiveness  of 
using  experimental  modal  testing  to  verify  the  adequacy  of  Finite  Element 
Models  (FEMs)  of  complex  missile  structures.  Before  pursuing  the  stated 
objective,  the  more  encompassing  question  of  "Why  FEMs?**  was  addressed. 

Insight  was  gained  by  examining  the  advantages  that  FEMs  offer  In  the  struc¬ 
tural  design  and  development  of  guided  missiles.  Basically,  missile  airframes 
and  supporting  hardware  are  designed  to  withstand  quasl-steady-state  loads 
(static  loads)  and  then  tested  and  modified  until  the  structures  withstand  the 
dynamic  loads  associated  with  real  world  deployment: 

A.  Static  Loads 


The  structural  designs  of  many  missiles  are  based  primarily  on  the 
quasi-steady  state  forces  or  the  use  of  static  accelerations,  g's,  for  the 
Inertia  loadings.  The  aerodynamic,  propulsive  (Internal  pressure  and  thrust), 
and  resulting  Inertia  loads  (propulsive  and  maneuvering)  are  considered  static 
forces.  Even  the  most  dynamic  of  loads  are  sometimes  approximated  with 
"equivalent  static  accelerations"  such  as  the  use  of  eight  g's  static  to  simu¬ 
late  the  loads  Imparted  to  a  missile  and  Its  launcher  when  a  helicopter 
crashes.  A  diagram  of  the  more  recent  development  of  the  FOG-M  structure  Is 
shown  In  Figure  1.  Even  this  simplified  or  truncated  design  approach  can  be 
enhanced  using  FEMs  early  In  the  development  process. 

Initially,  simple  beam  models  should  be  used  In  the  formulation  of  the 
mathematical  models  with  limited  use  of  plate  and  shell  elements.  Once  the 
models  are  established,  the  quasl-steady-state  or  static  forces  can  be  applied 
to  the  multinodes,  and  the  Internal  moment  and  shear  forces  derived  for 
detailed  stress  analysis.  If  more  sophisticated  depiction  of  Internal 
stresses  and  deflections  are  needed,  more  detailed  mathematical  models  can  be 
formulated  using  plate,  shell,  and  solid  elements  such  as  six  node  triangles 
and  eight  node  quadrilaterals.  Once  the  adequacy  of  the  simplified  model  has 
been  established  as  a  baseline,  the  extension  of  these  results  to  more 
sophisticated  models  Is  straightforward. 


B.  Dynamic  Forces 

The  use  of  FEMs  becomes  Imperative  If  d3mamlc  loads  are  used  In  the 
analysis.  The  basic  disciplines  are  well  established  to  Include  the  dynamic 
forces  and  the  related  structural  responses  early  In  the  design  process.  The 
prudent  approach  is  to  use  these  tools  which  will  significantly  reduce  costly 
redesigns,  engineering  change  proposals,  and  hardware  retrofits,  and  yield  more 
efficient  structures  earlier  in  the  development  sequence.  The  design  and 
development  process  will  remain  Iterative  as  prototype  hardware  becomes 
available  the  testing  of  actual  hardware  Is  possible.  The  difference  in  the 
structural  designs  is  based  on  actual  loads  and  not  in  oversimplified  static 
assumptions . 

Typical  missile  dynamic  loads  include  the  following: 

-shock  and  vibration  resulting  from  handling,  transportation,  and 
missile  launch  platform  dynamics  such  as  those  contributed  by 
helicopters  and  ground  vehicles. 

-launch  dynamics 

-impacts  associated  with  vehicle  Incidents  such  as  the  previously 
discussed  helicopter  crashes 

-blade  Induced  oscillations  from  helicopter  and  fixed  wing  aircraft 

-oscillatory  thrust  force 

-buffeting 

-acoustical  excitation 

-structural  response  resulting  from  thrusters 

-aerodynamic  force  resulting  from  missile  flexibility  or 
aeroelastlclty 

These  dynamic  loads  do  occur  in  real  life.  But,  in  the  "static  load” 
design  process,  the  stresses  associated  with  the  loads  are  hidden  in  the 
safety  factors  and/or  margins  of  safety  that  appear  in  the  stress  analysis. 

In  cases  where  these  margins  are  not  sufficient  to  cover  the  Ignored  dynamic 
loads,  the  faulty  designs  are  discovered  later  when  failures  occur  during 
testing  or  in  field  use.  In  many  cases,  the  dynamic  or  dynamic  related 
responses  are  actually  the  loads  that  ultimately  control  the  design  of  the 
missile.  Unfortunately,  the  loads  are  considered  later  in  the  development 
process  and  control  the  redesign  in  the  last  review.  As  shown  in  Figure  2, 
after  the  basic  development  process  is  complete,  flight  tests,  ground 
.ranspoi I .tlon  tests,  and  laboratory  tests  are  used  to  Identify  structural 
problems  resulting  from  dynamic  forcing  functions.  These  structural  problems 
are  corrected  and  an  iterative  process  follows  redesigning,  fabricating,  and 
testing  with  continual  modification  throughout  the  life  of  the  missile  system. 
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Figure  2.  Dynamic  force  test  approach. 


C.  Modal  Analysis  and  Testing 

The  need  for  adequate  mathematical  models  has  been  established  which 
substantiates  the  basic  argument  for  using  FEMs.  The  stated  purpose  of  this 
study  is  to  establish  that  good  engineering  correlation  can  be  obtained  bet¬ 
ween  the  FEM  predicted  modes  of  vibration  and  the  results  from  laboratory 
modal  tests.  Correlation  between  predicted  and  measured  modes  of  vibration 
establishes  the  basic  validity  of  the  models.  In  addition,  verification  of 
the  modes  proves  that  the  section  properties  and  mass  distributions  are 
correct.  Static  load  and  deflection  tests  ''an  and  should  be  used  in  the 
establishment  of  material  section  properties.  However,  the  modal  analysis  and 
testing  approach  carries  this  process  one  step  further  in  that  mass  distribu¬ 
tion  is  included  in  the  process.  Also,  the  importance  of  structural  features 
such  as  nonlinear  joints  and  structural  damping  are  an  inherent  part  of  the 
V  modal  analysis  process. 
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In  addition  to  verification  of  the  aatheutlcal  ■odela,  the  aodal  analysis 
and  testing  serve  a  significant  role  in  the  following  analyses: 

Free  Flight  Rockets  -identification  of  structural  damping  which 

controls  coupling  between  roll,  pitch,  and  yaw 

-aeroelastlclty  studies  for  the  long  L/D  low 
frequency  alaalles 

Guides  Missiles  -separation  of  structural  frequencies  froa 

guidance  and  control  loops 

-identification  of  ideal  bulkhead  location  for 
seeker  and  gyro  aountlng 

-Identification  of  actuator  locations 

D.  Missile  Under  Study 

A  prototype  aissile  structure  of  the  FOG-M  is  used  in  this  study*  In 
addition  to  providing  the  deaonstration  of  correlation  between  a  aatheaatlcal 
model  prediction  and  testing,  this  study  provides  a  coaparlson  between  the 
more  classical  normal  mode  testing  technique  and  the  now  more  popular  random 
vibration  approach.  Also,  observations  are  offered  concerning  the  strong 
points  and  shortcoalngs  of  the  two  test  techniques.  Finally,  Lessons  Learned 
in  the  performance  of  this  study  are  provided  for  future  consideration. 


II.  MISSILE  AIRFRAME  ANALYSIS 


A.  Body  Noraal  Mode  Measureaentt 

Two  accelerometers  were  first  placed  at  posit  loos  I  and  2  shotra  In 
Figure  3  with  the  exciter  coupled  to  the  bottom  of  the  missile  directly  below 
accelerometer  position  2.  Results  are  shown  In  Figure  A  with  the  major  reso¬ 
nances  being  283.6  Hz  at  location  I  and  132.6  Hz  at  location  2.  The  axes  of 
these  plots  are  linear  peak  g's  versus  frequencies  from  15  to  500  Hz.  Also 
note,  the  minor  resonances  that  fall  between  85  and  95  Hz  for  both  locations 
which  will  be  shown  later  as  resonances  associated  with  other  subsystems 
Internal  to  the  missile. 

Figure  5  shows  a  repeat  of  the  first  sine  sweep  (Figure  4  data)  which 
shows  acceptable  repeatability.  The  cross-hair  must  be  aligned  with  one  of 
the  resonances.  This  alignment  was  arbitrarily  selected  as  131.4  Hz  for 
accelerometer  2.  As  shown,  the  cross-hair  Indicates  a  shift  for  the  same  reso¬ 
nance  at  accelerometer  location  1  (bottom  plot).  This  shift  Is  within  the 
resolution  of  this  Instrumentation  system.  The  shift  represents  a  difference 
of  3  to  4  percent  from  the  indicated  frequency  of  131.4  Hz. 

Next,  the  transducer  located  at  position  1  was  moved  to  location  3 
and  the  sine  sweep  repeated.  The  driving  acceleration  was  somewhat  less  as 
shown  by  comparing  location  2  peak  g's  on  Figure  4  and  Figure  6  but  is  within 
acceptable  tolerances.  The  lower  frequency,  128.4  Hz,  (same  resonance  as  the 
131.4  Hz  previously  noted)  was  again  evident  at  locations  2  and  3  and  a  second 
resonance  of  289.7  Hz  was  prevalent  at  location  3  which  Is  close  to  the  283.6 
Hz  observed  during  the  first  test. 

The  last  sine  sweep  test  conducted  was  with  a  transducer  at  location 
2  and  a  moveable  accelerometer  at  location  4  on  the  Inside  lip  of  the  can¬ 
tilevered  fiber  optics  spool.  The  results  are  shown  In  Figure  7  with  the  only 
predominant  resonance  at  location  4  being  at  approximately  92  Hz.  Additional 
measurements  at  various  locations  along  the  length  of  the  spool  would  have 
been  helpful  but  would  have  necessitated  disassembly  of  the  structure.  Other 
measurements  were  made  on  the  outer  skin  at  the  approximate  location  of  the 
base  of  the  spool,  and  accelerations  were  small.  Consequently,  It  is  safe  to 
Infer  that  the  first  cantilever  mode  for  the  spool  is  92  Hz. 

Now,  recall  the  resonance  observed  on  the  front  portion  of  the 
missile  at  a  frequency  between  85  and  95  Hz.  An  earlier  modal  experiment  was 
conducted  on  the  TV  seeker  and  glmbal  mounting  and  Is  documented  In  Reference 
1.  The  first  resonance  of  the  cantilever  TV  seeker  was  found  to  be  between  95 
and  105  Hz.  The  exact  glmbal  system  and  seeker  were  probably  a  later  version 
than  the  system  examined  In  this  study,  but  the  range  of  the  frequencies 
strongly  suggests  that  both  resonances  are  cantilevered  beam  modes  of  the 
mounted  seeker.  During  these  more  recent  tests,  the  missile  was  not 
disassembled,  but  It  appears  conservative  to  assume  that  the  first  cantilever 
mode  of  the  seeker  Is  approximately  90  Hz. 
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The  next  step  In  the  oormel  aode  aeeeureaent  scheae  vas  to  conduct 
dwell  tests  at  Che  frequencies  Identified  In  the  sine  sweep  tests.  The 
exciter  was  placed  In  exactly  the  sase  poaitlon  as  for  the  sine  sweep  test. 

One  of  the  two  available  acceleroaeters  was  specified  as  the  fixed  stationary 
reference  to  ensure  the  saae  relative  acceleroaeter  la  being  applied  for  each 
test  setup  and  to  establish  phase  relationships  for  the  aeasureaents.  In 
these  measurements,  the  reference  accelerometers  were  place  above  the  excita¬ 
tion  source  or  location  5  in  Figure  8.  Also  shown  on  this  figure  are  the  nine 
locations  for  the  roving  accelerometers.  The  driving  force  Is  then  set  at  a 
frequency  near  the  resonances  Identified  earlier  and  Halted  frequency  sweeps 
are  run  up  and  down  to  pinpoint  the  resonance  to  be  studied.  The  first  reso¬ 
nance  was  between  128.4  Hz  and  132.6  Hz  from  the  sine  sweep  testa  depending  on 
location  of  the  measurement.  In  the  case  of  these  limited  frequency  sweeps, 
the  resonant  frequency  was  identified  as  130.3  Hz  or  between  the  previous 
measurements.  Once  this  measurement  scheme  was  established,  data  were 
recorded  In  a  format  similar  to  the  first  four  columns  of  Table  1. 


TABLE  1.  Airframe  First  Body  Bending 
FREQOENCY  -  130.3  Hz 


LOC 

AMP 

REF  AMP 

PHASE 

C0S(»i-d2) 

SHAPE 

1 

.507 

.738 

-158.9 

+  .99 

+  .50 

2 

.221 

.753 

-140.0 

+  .98 

+  .22 

3 

.349 

.751 

4.6 

-  .92 

-  .32 

4 

.659 

.748 

3.6 

-  .91 

-  .60 

5 

.688 

.738 

2.8 

-  .91 

-  .63 

6 

.595 

.729 

-  5.5 

-  .84 

-  .50 

7 

.517 

.747 

-  16.9 

-  .71 

-  .37 

8 

.342 

.753 

-  81.7 

+  .33 

+  .11 

9 

1.390 

.794 

0 

+1,00 

+1.00 

Next,  the  maximum  acceleration  was  identified  which  was  at  location  9 
for  this  particular  frequency.  Then,  the  cosine  of  the  angle  between  the 
maximum  acceleration  and  the  phase  angle  at  each  of  the  ocher  eight  locations 
was  computed  with  the  results  shown  in  column  5  of  the  table.  Next,  the  eight 
remaining  acceleratluu  measurements  were  normalized  to  the  maximum  accelera¬ 
tion  and  multiplied  by  the  cosine  of  the  phase  angle  differences.  This  pro¬ 
duct  becomes  the  mode  shape.  This  final  mode  shape  will  be  shown  later  after 
a  discussion  of  the  NASTRAN  modal  predictions. 


This  same  procedure  was  then  repeated  for  the  second  noted  resonance 
which  was  around  283.6  Hz  to  289.7  Hz  from  the  sine  sweep  tests.  The  limited 
frequency  sweep  yielded  a  resonance  of  280.0  Hz  again  between  the  earlier 
measurements.  The  sine  dwell  tests  yield  the  data  In  Table  2  which  again  will 
be  shown  In  graphical  format  after  the  following  discussion  on  the  NASTRAN 
modal  predictions. 


TABLE 

2.  Airframe  Second  Body  Bending 
FREQUENCY  -  280.0  Hz 

LOG 

AMP 

REG  AMP  PHASE 

C0S((>T-»2) 

SHAPE 

1 

.697 

.404 

-133.0 

+  .74 

+  .24 

2 

.376 

.414 

-110.0 

+  .41 

+  .07 

3 

.326 

.414 

-  53.0 

-  .54 

-  .08 

4 

.361 

.417 

-  16.6 

-  .93 

-  .16 

5 

.418 

.416 

-  .5 

-1.0 

-  .20 

6 

.419 

.417 

+  1.0 

-1.0 

-  .19 

7 

.381 

.417 

-  3.2 

-  .99 

-  .18 

8 

.153 

.419 

-  69.8 

-  .27 

-  .02 

9 

2.110 

.416 

-175.5 

+1.0 

+1.00 

* 

.375 

-169.0 

** 

.495 

-155.0 

right  side  vertical  direction 

location  9 

**  bottom  side  vertical  direction  location  9 


B.  NASTRAN  Modal  Analysis 

A  beam  model  for  the  FOG-H  Is  shown  In  Figure  9.  When  three  dimen¬ 
sional  modal  plots  are  obtained  from  the  NASTRAN  plotting  routines.  It  is 
extremely  difficult  to  ascertain  the  higher  frequency  body  bending  modes 
because  the  relatively  low  frequency  wings  (20  Hz  range)  and  fins  bending 
modes  (45  Hz  range)  are  extrem^y  flexible  and  the  displacements  are  nor¬ 
malized  to  the  wings  and  fins  even  though  the  mode  of  vibration  Is  primary 
body  bending.  An  Illustration  of  this  point  Is  shown  In  Figure  10  which  Is  a 
plot  of  the  NASTRAN  16th  mQde  In  only  the  x-z  plane.  Even  though  the  mode  is 
primarily  body  bending,  the  wing  and  fin  displacement  overshadow  the  displace¬ 
ments  of  Interest.  This  observation  Is  even  more  prevalent  when  trying  to 
interpret  three  dimensional  NASTRAN  mode  plots. 
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In  order  to  resolve  the  NASTRAN  predicted  bending  modes  for  the  body 
bending  from  those  of  the  wings  and  fins,  the  mathematical  model  was  then 
recomputed  with  the  wings  and  fins  totally  removed  and  then  with  lumped  masses 
In  place  of  the  appendages  at  the  Junctures  with  the  missile  body.  Also,  hand 
calculations  were  made  using  average  section  properties  and  assuming  equally 
distributed  mass  for  a  simple  free-free  beam  for  the  case  when  the  wings  and 
fins  are  removed.  Results  are  shown  In  Figure  11  and  the  following  table: 


TABLE  3.  Calculated  First  Body  Bending 


CONDITION 


FREQUENCY 


First  body  bending  167.3  Hz 

w/o  wings  &  fins 
maniial  calculations 

First  body  bending  160.8  Hz 

w/o  wings  &  fins 

NASTRAN 

First  body  bending 

Itimped  wing  &  fin  masses  14A.7  Hz 

NASTRAN 

First  body  bending  147.2  Hz 

distributed  wing  &  fin  masses 

NASTRAN 

Second  body  bending  467.77  Hz 

distributed  wing  &  fin  masses 


This  exercise  provided  confidence  that  the  NASTRAN  mathematical  model 
is  approximately  correct  and  the  modes  have  been  correctly  Identified  when 
masked  by  the  large  displacement  associated  with  the  lower  frequency  wing  and 
fin  modes  of  vibration.  Now  that  the  background  has  been  provided  on  the 
mathematical  model,  the  comparison  of  NASTRAN  modal  predictions  and  the 
measurements  made  using  the  normal  modal  analysis  scheme  can  continue. 
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C.  Comparison  Measureaenta  .  id  NASTRAN  Modal  Analysis 

Figure  12  shows  comparison  of  the  predicted  and  measured  first  mode  of 
vibration.  The  NASTRAN  data  was  extracted  from  the  computer  printout,  only 
the  body  bending  was  plotted  and  the  data  renormalized  to  the  maximum  peak  for 
the  body  elements.  Excellent  correlation  was  obtained  for  the  mode  shape 
which  leaves  little  doubt  that  the  correct  modes  are  being  compared. 

Obviously,  a  discrepancy  exists  between  the  predicted  and  measured 
frequencies,  l.e.,  147  Hz  predicted  and  130  Hz  measured.  One  possible  source 
of  error  In  the  mathematical  model  Is  the  numerous  structural  joints  In  the 
missile.  There  are  seven  structural  Joints  In  the  FOG-M  missile  body.  Each 
of  these  joints  can  Introduce  flexibility  not  In  the  basic  shell  from  the 
following: 

-bolt(s)  bending 

>def lection  In  shell  due  to  beam  bending 

-beam  deflection 

-deflection  due  to  clearances  (slop) 

A  possible  refinement  to  the  mathematical  model  In  missile  body  simula¬ 
tion  would  be  to  add  additional  elements  for  the  joints  and  study  the  effect  on 
modal  predictions.  This  refinement  or  others  are  necessary  If  the  difference 
between  130  and  147  Hz  la  unacceptable. 

The  second  measured  resonance  was  260  Hz.  A  plot  of  the  measured  mode 
shape  and  the  NASTRAN  second  mode  of  vibration  Is  shown  In  Figure  13.  This  mode 
is  definitely  not  the  actual  second  mode  and  occurs  at  a  frequency  significantly 
below  what  would  be  expected  based  on  a  first  mode  of  130  Hz.  The  second  mode 
should  be  at  a  frequency  which  is  approximately  two  and  three  quarters  to  three 
times  the  lower  frequency  or  358  to  390  Hz.  This  observation  Is  further 
substantiated  by  the  fact  that  the  vertical  vibration  was  also  measured  on  one 
side  and  on  the  bottom  of  the  missile  as  shown  in  Table  2.  The  magnitude  was 
significantly  less.  Indicating  that  missile  body  bending  was  not  taking  place. 
This  resonance  Is  attributed  to  some  localized  deflections  and  Is  not  of  signi¬ 
ficant  Importance  In  understanding  the  overall  response  of  the  missile. 
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III.  WING(S)  ANALYSIS 


A.  Left  Wing  Normal  Mode  Measurements 

The  FOG-M  has  four  wings.  The  configuration  chosen  for  this  modal 
test  was  to  position  one  set  of  wings  horlzonally.  The  left  wing  refers  to 
the  left-hand  side  when  facing  the  direction  of  the  forward  motion  of  the 
missile.  The  configuration  tested  actually  had  three  wings  since  the  bottom 
wing  was  removed  so  the  exciter  could  be  positioned  closer  to  the  center  line 
of  the  wing.  As  will  be  shown  later,  all  three  wings  could  have  been  measured 
slmxiltaneotisly  at  each  identified  modal  frequency.  This  sequence  would  have 
been  directly  comparable  to  that  used  for  the  NASTRAN  modal  predictions. 
However,  available  instrumentation  and  time  suggest  another  approach.  The 
same  information  can  be  gained  from  the  less  time  consuming  task  of  measuring 
one  wing.  This  is  the  approach  followed  in  this  series  of  tests  except  that 
in  addition  to  the  left  wing  which  is  free  of  voids,  cutouts,  etc.,  the  right 
wing  was  measured  separately  to  study  the  effect  of  a  cutout  in  the  top  of  the 
wing.  The  cutout  is  used  to  show  internal  moving  wing  parts  such  as  springs 
and  rods  in  the  mock-up  missile.  A  schematic  of  the  top  view  of  the  missile 
is  shown  in  Figure  14.  The  electrodynamic  shaker  is  also  shown  in  the  figure 
placed  directly  beneath  the  missile  and  connected  to  the  underside  of  the 
missile  with  a  short  rod  or  stinger.  A  tapped  hole  in  a  bulkhead  of  Che 
missile  body  conveniently  accepted  the  stinger. 

The  sequence  was  to  again  conduct  sine  sweeps  over  the  frequency 
bands  of  interest.  The  results  from  these  sine  sweeps  are  shown  In  Figure  15 
(15  to  200  Hz),  Figure  16  (150  to  400  Hz),  and  Figure  17  (350  to  500  Hz).  The 
Cop  curve  In  each  figure  Is  the  accelerometer  located  on  the  tip  of  the  wing 
(position  1)  and  the  bottom  plot  on  each  figure  Is  the  accelerometer  located 
at  the  base  of  the  wing  (position  6).  Examination  of  these  plots  shows  the 
major  resonances  to  be  20  Hz,  173  Hz,  191  Hz,  and  506  Hz  (not  shown).  Each 
one  of  the  resonances  was  then  Investigated  in  detail  as  presented  In  the 
following  paragraphs. 

The  resonances  were  again  Isolated  by  conducting  limited  frequency 
band  sweeps  near  Che  Indicated  frequencies  until  the  peak  response  was 
Isolated.  The  first  resonance,  19.6  Hz,  data  are  given  In  Table  4,  and  the 
mode  plot  provided  In  Figure  18  with  the  predicted  NASTRAN  modal  frequency  and 
shape.  Almost  perfect  correlation  was  obtained  between  Che  NASTRAN  mode  shape 
and  Che  measured  modal  response. 


TABLE  4. 

Left  Wing 

First 

Bending  Mode 

FREQUENCY 

-  19.6 

Hz . 

LOG 

AMP 

REF  AMP 

PHASE 

C0S(di-d2) 

SHAPE 

Ktlp) 

2.90 

.29  -144.6 

+1.00 

+1.00 

2 

2.00 

.28 

-141.0 

+1.00 

+  .69 

3 

1.30 

.29 

-139.0 

+1.00 

+  .45 

4 

.83 

.28 

-131.0 

+  .97 

+  .28 

5 

.36 

.31 

-112.0 

+  .84 

+  .10 

6(base 

)  .22 

.28 

-  25.0 

-  .49 

-  .04 
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The  next  two  modes  measured  were  second  mode  bending  modes  of  175 
and  192  Hz.  The  measured  mode  shapes  are  given  In  Tables  5  and  6.  A  graphi¬ 
cal  comparison  of  the  two  modes  is  given  in  Figure  19  which  shows  both  modes 
to  be  second  mode  beam  bending.  Next,  these  two  measured  modes  are  plotted 
versus  the  NASTRAN  mode  shapes  in  Figures  20  and  21. 


TABLE  5.  Left  Wing  Second  Bending  Mode 
FREQUENCY  -  175  Hz 


LOG 

AMP 

REF  AMP 

PHASE 

C0S((ti-<»?) 

SHAPE 

Ktlp) 

9.40 

.29 

65.0 

-  .96 

-  .90 

2 

2.70 

.28 

-71.0 

+  .48 

+  .13 

3 

9.80 

.34 

233.0 

+1.00 

+  .98 

4 

10.00 

.38 

228.0 

+1.00 

+1.00 

5 

5.60 

.30 

249.0 

+  .93 

+  .52 

6 (base) 

.39 

.28 

-41.0 

-  .02 

0 

TABIE 

6.  Left  Wing  Second  Bending 

FREQUENCY  -  192 

Hz 

LOG 

AMP 

REF  AMP 

PHASE 

C0S(^l-«2) 

SHAPE 

Ktlp) 

7.30 

.41 

33.5 

+1.00 

+1.00 

2 

1.80 

.30 

-  93.0 

-  .59 

-  .15 

3 

6.30 

.40 

215.0 

-1.00 

-  .86 

4 

6.50 

.43 

209.0 

-1.00 

-  .89 

5 

4.60 

.31 

237.0 

-  .92 

-  .58 

6 (base) 

.29 

.28 

-  15.0 

+  .66 

+  .03 

Again,  good  correlation  was  obtained  as  to  mode  shapes.  However, 
the  frequencies  do  not  agree.  The  theoretical  frequencies  were  90  Hz,  119  Hz, 
119  Hz,  and  120  Hz.  See  Figure  22  for  predicted  mode  shapes.  There  are  four 
frequencies  because  of  the  four  wings  which  provide  four  coupled  eigenvectors 
and  eigenvalues  for  each  mode.  Again,  the  two  measured  frequencies  for  this 
one  wing  were  175  and  192  Hz. 

The  first  candidate  source  that  could  explain  the  difference  bet¬ 
ween  measured  and  predicted  higher  bending  mode  frequencies  is  the  fixity  at 
the  base  of  the  wings  and  the  flexibility  of  the  missile  body  In  the  region  of 
the  wlng(s)  attachment.  The  model  could  be  modified  to  include  torsional 
springs  and/or  simulants  for  nonlinear  effects  in  an  attempt  to  raise  the 
NASTRAN  higher  mode  frequencies  while  not  significantly  modifying  the  first 
bending  mode.  This  provided  excellent  correlation  for  both  frequency  and 
mode  shape.  Another  possible  explanation  for  the  higher  frequency  difference 
could  be  the  mass  distribution  used  in  the  models  of  the  wings.  The  mass  of 
the  internal  working  parts  of  the  wings  was  distributed  along  the  wing  which 
would  tend  to  lower  the  predicted  natural  frequencies. 
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The  additional  modeling  refinements'  are  not  warranted  for  the 
majority  of  missile  analyses.  The  displacements  which  control 
strains /stresses  are  a  function  of  frequency  squared.  The  ratio  of  the 
squares  of  the  first  and  second  measured  modes  for  the  wings  Is  19.62/1752  or 
1.25  percent.  The  Power  Spectral  Density  (PSD)  levels  for  the  vast  majority 
of  the  vibration  environments  (forcing  functions)  actually  decrease  In  magnl- 
tude  (g2/Hz)  at  the  higher  frequencies  (see  Reference  2).  This  observation  Is 
readily  understood  when  the  natural  frequency  of  the  carrier  vehicle  and 
mounting  structure  Is  taken  Into  consideration.  Typically,  the  suspension 
system  will  be  In  the  lower  frequency  range;  consequently,  a  reduction  In 
vibration  energy  at  the  higher  frequencies  should  be  expected.  Also,  the 
transportation  and  missile  delivery  platform  vibrations  are  more  severe  than 
flight  vibrations  such  as  acoustic  and  buffeting.  The  ground  based  environ¬ 
ments  control  the  structural  design  for  vibrations. 

The  last  measured  wing  mode  was  located  at  506  Hz  and  the  measured 
values  are  given  In  Table  7.  The  measured  mode  shape  and  NASTRAN  prediction 
are  given  In  Figure  23.  Just  as  with  the  second  bending  mode,  good  mode  shape 
correlation  was  obtained,  but  a  wider  dispersion  between  measured  and  pre¬ 
dicted  frequencies  was  obtained,  e.g.,  506  Hz  measured  with  the  four  predicted 
frequencies  falling  between  278  and  319  Hz.  Since  the  second  bending  mode 
frequency  did  not  agree,  even  a  leaser  degree  of  agreement  Is  expected  when 
comparing  frequencies  for  the  third  bending  mode.  The  previously  mentioned 
differences  In  fixity  and  mass  distributions  could  explain  the  difference  In 
frequencies . 


TABLE 

7.  Left  Wing  Third  Bending  Mode 

FREQUENCY  -  506 

Hz 

LOG 

AMP 

REF  AMP 

PHASE 

C0S(»i-<>7) 

SHAPE 

l(tlp) 

6.30 

.46 

-120.0 

+  .99 

+  .52 

2 

7.70 

.54 

16.0 

-  .82 

-  .53 

3 

.79 

.50 

-240.0 

-  .36 

-  .02 

4 

6.30 

.50 

-110. 0 

+  .95 

+  .50 

5 

11.90 

.40 

-129.0 

+1.00 

+1.00 

6 (base) 

.77 

.54 

-  6.6 

-  .54 

-  .03 
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B.  Left  Wing  Random  Mode  Measurements. 

The  procedure  for  conducting  the  random  modal  measurements  Is  given 
In  Reference  3.  In  summary,  the  procedure  calls  for  exciting  the  structure 
with  random  white  noise  (flat  spectrum)  energy  and  recording  the  measurements 
for  use  In  subsequent  analyses.  The  preferred  technique  to  Isolate  which  fre¬ 
quencies  to  study  Is  to  measure  the  ratio  of  the  accelerometer  of  Interest  and 
the  Input  force  In  units  of  g's/lbf.  Plots  for  the  six  positions  defined  In 
Figure  14  are  given  In  Figures  24  through  29.  In  this  format,  the  Iden¬ 
tification  of  which  natural  frequencies  to  study  Is  difficult;  consequently, 
the  equipment  manufacturer  provides  a  code  called  "find  mode"  which  provides  a 
weighted  average  of  the  g's/lbf  for  all  locations  under  consideration.  This 
presentation  for  the  frequency  band  from  15  to  210  Hz  Is  given  In  Figure  30. 
The  predominant  frequencies  for  this  band  are  21,  173,  and  195  Hz  which  pro¬ 
vide  excellent  agreement  with  values  measured  using  the  normal  mode  method. 

In  addition,  the  first  and  second  beam  bending  mode  shapes  show  acceptable 
agreement  with  the  previous  modal  survey  technique  (see  Figures  31,  32,  and 
33).  Unfortunately,  other  frequencies  shown  In  Figure  30  might  also  warrant 
study.  For  example,  the  fourth  highest  amplitude  Is  approximately  119  Hz. 

The  software  has  no  method  to  differentiate  good  modes  from  nonmodes  of  other 
parts  of  the  missile  resonating  such  as  Internal  missile  components.  This 
lack  of  ability  to  differentiate  good  from  superfluous  modes  Is  amplified  when 
the  next  step  Is  taken  In  the  software  to  measure  the  mode  shape.  Figure  34 
shows  a  mode  plot  of  the  119-120  Hz  frequency.  .  The  mode  appears  to  be  that  of 
a  plnned-plnned  beam  and  not  of  the  cantilevered  wing.  Assuming  that  the 
signal  Is  out  of  the  noise,  two  sources  could  explain  this  apparent  mode. 

First ,  the  measurement  scheme  was  designed  to  measure  beam  bending  with  the 
accelerometers  located  along  the  center  line  of  the  symmetrical  wing.  The 
mode  could  be  an  artifact  of  the  first  torsional  mode  which  should  fall  be¬ 
tween  100  and  200  Hz  with  the  first  beam  bending  mode  of  20  Hz.  Placement  of 
accelerometers  in  a  matrix  along  the  plane  of  the  wing  would  have  allowed  for 
the  Identification  of  torsional  modes.  This  modeling/ test  program  was 
Initiated  as  a  simple  NASTRAN  beam  model  approach;  consequently,  modeling  and 
measuring  the  torsional  modes  of  vibration  were  out  of  scope  of  the  project. 
Another  explanation  of  the  120  Hz  mode  Is  that  some  other  component,  missile 
substructure,  appendage,  etc..  Is  vibrating  and  the  resulting  wing  mode  shape 
falls  between  the  first  and  second  beam  bending  modes.  In  the  future,  great 
care  should  be  taken  to  ensure  that  actual  modes  are  being  measured  and  not 
artifacts  of  other  natural  resonances  when  using  the  Random  Mode  Measurements 
technique.  The  slower  Normal  Mode  Measurements  method  requires  that  each 
apparent  mode  of  vibration  be  Identified  and  verified  as  a  predominate  struc¬ 
tural  mode  as  the  tester /analyst  moves  up  the  frequency  scale. 

C.  Right  Wing  Normal  Mode  Measurements. 

The  right  wing  normal  modes  were  measured  by  the  same  procedure  pre¬ 
viously  described  for  the  left  wing.  The  locations  of  accelerometers  and  the 
electrodynamic  shaker  are  given  in  Figure  35.  The  primary  differences  between 
the  right  wing  and  left  wing  are  the  presence  of  the  previously  described 
cutout  and  a  longer  rod  or  stinger  to  attach  the  shaker  to  the  missile.  The 
sine  sweep  data  are  given  In  Figure  36  (15  Hz  to  200  Hz),  Figure  37  (150  Hz  to 
400  Hz),  and  Figure  38  (350  Hz  to  600  Hz).  The  identified  resonant  frequen¬ 
cies  are  21  Hz,  160  Hz,  433  Hz,  and  481  Hz. 
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As  with  the  other  missile  substructures,  sine  dwell  tests  were  con¬ 
ducted  at  these  frequencies.  The  data  are  given  In  Table  8  for  the  21  Hz  mode, 
and  the  cantilevered  beam  mode  shape  is  provided  in  Figure  39.  Again, 
excellent  correlation  Is  provided  for  the  left  wing  first  mode  both  with  data 
from  the  normal  mode  method  and  the  random  mode  method,  and  with  the  NASTRAN 
model  mode  predictions . 


TABLE  8.  Right  Wing  First  Bending  Mode 
FREC^JEMCY  -  21.0  Hz 


LOG 

AMP 

REF  AMP 

PHASE 

C0S((tT-<>?) 

SHAPE 

l(tip) 

2.50 

.26 

-167.0 

+1.00 

+1.00 

2 

1.62 

.23 

-167.0 

+1.00 

+  .65 

3 

1.08 

.22 

-166.0 

+1.00 

+  .43 

4 

.52 

.21 

-164.0 

+1.00 

+  .21 

5 

.13 

.20 

-140.0 

+  .89 

+  .05 

6 (base) 

.17 

.21 

-  9.0 

-  .93 

-  .06 

The  second  indicated  mode  from  the  sweep  tests  was  a  strong  mode  at 
160  Hz.  The  dwell  tabulated  data  are  given  in  Table  9  and  plots  in  Figure  40. 
This  Is  definitely  one  of  the  second  beam  bending  modes.  The  frequency  Is 
lower  than  the  lowest  measured  comparable  mode  for  the  left  wing  of  175  Hz. 

Of  course,  both  the  160  and  175  Hz  measured  frequencies  are  higher  than  the 
NASTRAN  model  prediction  of  90  Hz.  The  differences  between  the  frequencies  of 
right  and  left  wing  second  bending  modes  are  probably  due  to  the  effect  of  the 
cutout.  If  this  is  true,  some  calibration  Is  provided  as  to  what  could  be 
expected  when  using  modal  testing  to  study  composite  structure  defects.  The 
mode  shapes  for  the  right  and  left  wings  are  similar,  but  the  frequencies  are 

different,  l.e.,  160  Hz  versus  175  Hz.  The  effect  of  the  cutout  is 

detectable;  however,  the  relative  size  of  the  defect  (cutout)  is  large  com¬ 
pared  to  the  wing  plan  form  area.  These  observations  provide  some  insight 

Into  using  modal  testing  to  detect  flaws  in  composite  structures.  Future 

Investigators  should  give  attention  to  the  second  and  higher  modes  of  vibra¬ 
tion  when  studying  Internal  flaws. 


TABLE  9.  Right  Wing  Second  Bending  Mode 
FREQUENCY  -  160  Hz 


LOG 

AMP 

REF  AMP 

PHASE 

cos(di~*^7 ) 

SHAPE 

Ktlp) 

6.11 

.35 

+  41.0 

+1.00 

-1.00 

2 

2.01 

.32 

+243.0 

-  .93 

+  .31 

3 

5.90 

.35 

+228.0 

-  .99 

+  .96 

4 

4.75 

.36 

+208.0 

-  ,97 

+  .75 

5 

3.95 

.30 

+250.0 

-  .87 

+  .56 

6 (base) 

.27 

.31 

-  21.0 

+  .47 

-  .02 
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Minor  responses  were  observed  at  172  and  214  Hz  In  the  sine  sweep 
data.  The  relative  magnitudes  are  Insufficient  to  Indicate  major  bending 
modes  of  vibration. 

The  last  two  modes  to  be  studied  on  the  right  wing  were  modes  located 
at  433  Hz  and  481  Hz.  Tabulated  data  are  given  in  Tables  10  and  11,  and  the 
mode  plots  are  provided  In  Figures  41  and  42.  The  modes  appear  to  be 
approaching  the  third  beam  mode  measured  on  the  left  wing  which  was  a  classi¬ 
cal  third  bending  mode  at  506  Hz.  The  poor  definition  of  the  right  wing  mode 
shape  Is  probably  again  due  to  the  cutout.  The  defect  evidently  has  more 
Influence  at  the  higher  modes. 


TABLE  10.  Right  Wing  Bending  Mode 
FREQUENCY  -  433  Hz 


LOG 

AMP 

REF  AMP 

PHASE 

C0S(>>T()>2) 

SHAPE 

l(tlp) 

5.10 

2.20 

+  90.0 

-  .92 

-  .47 

2 

10.00 

1.20 

-  67.0 

+1.00 

+1.00 

3  • 

4.70 

1.40 

-  6.0 

+  .48 

+  .23 

4 

3.50 

1.40 

+  5.0 

+  .31 

+  .11 

5 

6.50 

1.20 

+  19.0 

+  .70 

+  .46 

TABLE  11.  Right  Wing  Bending  Mode 
FREQUENCY  -  481  Hz 


LOC 

AMP 

REF  AMP 

PHASE 

C0S(»-i-»2) 

SHAPE 

Ktlp) 

5.80 

.71 

+160.0 

-  .87 

-  .34 

2 

14.80 

.95 

+  9.8 

+1.00 

+1.00 

3 

5.60 

.89 

+158.0 

-  .85 

-  .32 

4 

2.58 

.88 

+267.0 

-  .22 

-  .04 

5 

2.37 

.72 

+239.0 

-  .65 

-  .io 

6 (base) 

.66 

1.00 

+  72.0 

+  .47 

+  .02 

IV.  FIN(S)  ANALYSIS 


The  four  fins  were  connected  to  the  missile  structure  via  linkages,  pins, 
actuators,  and  other  hardware.  These  structures  contained  a  significant  amount 
of  clearance  (slop).  During  flight,  these  mechanisms  were  energized  with 
pneumatic  pressure  and  the  greater  portion  of  these  nonlinearities  were 
eliminated.  One  approach  to  these  tests  would  have  been  to  apply  the  appli¬ 
cable  pressures  and  to  pin  or  fix  any  remaining  joints  that  might  introduce 
excessive  clearance.  During  these  particular  tests,  a  simplified  approach  was 
followed  where  two  fins  were  connected  to  either  side  of  a  .625  Inch  cube 
block.  The  cube  block  was  connected  to  the  exciter  with  a  stinger  (small 
diameter  threaded  rod)  and  a  force  transducer.  A  schematic  of  the  test  setup 
Is  shown  In  Figure  43.  The  length  of  the  rod  connecting  the  fins  to  the  cube 
block  was  .191  Inch  which  Is  representative  of  the  length  In  the  actual 
missile  application  and  Is  the  same  length  chosen  In  the  NASTRAN  math  model. 
Again,  normal  mode  measurements  and  random  response  techniques  were  used  to 
Identify  the  modes  of  vibrations  for  comparison  to  NASTRAN  predictions. 

A.  Fin  Normal  Mode  Y-2  Direction. 

The  first  normal  mode  measurements  were  made  In  the  Y-Z  directions. 

A  plot  of  the  sine  sweep  data  Is  shown  In  Figure  44.  The  vertical  axis 
depicts  transmlsslblllty,  that  Is,  the  acceleration  at  the  right  hand  tip 
(accelerometer  1)  divided  by  the  acceleration  at  the  top  of  the  cube  mounting 
block.  As  given  In  the  figure,  the  highest  amplitude  resonances  were  42,  291, 
320,  and  470  Hz.  Next,  each  of  these  resonances  were  studied  In  detail. 

Limited  band  frequency  sweeps  were  conducted  at  a  nominal  42  Hz  frequency  and 
the  peak  response  frequency  was  established  to  be  43  Hz.  Dwell  tests  were 
then  conducted  and  the  results  are  given  In  Table  12  and  Figure  45. 


TABLE  12.  Fin  First  Bending  Mode 
FREQUENCY  -  43  Hz 


LOG 

AMP 

REF  AMP 

PHASE 

COS  (01  •*0?  ) 

SHAPE 

l(tip) 

10.37 

.84 

-127.0 

+1.00 

+1.00 

2 

8.40 

.60 

-  92.0 

+  .81 

+  .66 

3 

5.99 

.61 

-  67.0 

+  .50 

+  .29 

4 

4.04 

.63 

-  57.0 

+  .34 

+  .13 

5 

1.97 

.64 

-  52.0 

+  .25 

+  .05 

6(base) 

.71 

.64 

-  14.0 

-  .39 

-  .03 

A 


Also  shown  on  Figure  45  Is  the  NASTRAN  predicted  mode  of  43 
Hz.  Both  the  frequency  and  modes  shape  correlations  are  considered 
acceptable.  The  model  could  be  changed  to  provide  exact  correlation 
if  desired.  The  next  three  resonances  studied  were  294,  320,  and  328 
Hz  which  bracket  the  two  second  bending  modes  of  the  fins.  Tables  13, 
14,  and  15  sunnarlze  these  data. 
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TABLE 

13.  Fin 

Second  Bending  Mode 

FREQUENCY  - 

294  Hz 

LOC 

AMP 

REF  AMP 

PHASE 

C0S(»i-»?) 

SHAPE 

l(tlp) 

10.54 

1.94 

+126.2 

+1.00 

+1.00 

2 

2.09 

2.07 

+144.9 

+  .95 

+  .19 

3 

2.36 

2.07 

«  8.3 

-  .70 

-  .16 

4 

5.33 

2.02 

-  14.3 

-  .77 

-  .39 

5 

4.92 

2.03 

-  12.1 

-  .75 

-  .35 

6 (base) 

3.07 

2.05 

-  5.1 

-  .66 

-  .19 

TABLE 

14.  Fin  Second  Bending  Mode 

FREQUENCY  - 

320  Hz 

LOC 

AMP 

REF  AMP 

PHASE 

C0S((>i-«2) 

SHAPE 

l(tlp) 

13.78 

2.22 

+  72.7 

+1.00 

+1.00 

2 

.95 

2.13 

+117.5 

+  .71 

+  .05 

3 

5.47 

1.93 

-  32.1 

-  .26 

-  .10 

4 

8.89 

2.00 

-  50.3 

-  .54 

-  .35 

5 

8.50 

1.90 

-  35.7 

-  .32 

-  .20 

6 (base) 

3.64 

1.95 

-16.90 

+  .01 

+  .00 

TABLE 

15.  Fin 

Second  Bending  Mode 

FREQUENCY  - 

328  Hz 

LOC 

AMP 

REF  AMP 

PHASE 

C0S(<>i-'*2) 

SHAPE 

l(tlp) 

13.50 

2.30 

+  62.5 

+1.00 

+1.00 

2 

3.79 

1.94 

+105.8 

+  .72 

+  .20 

3 

5.22 

1.87 

-  63.9 

-  .56 

-  .22 

4 

10.01 

1.94 

-  78.5 

-  .78 

-  .58 

5 

7.51 

1.79 

-  51.2 

-  .40 

-  .22 

6 (base) 

3.07 

1.81 

-  24.3 

+  .06 

+  .01 

The  mode  shapes  are  shown  In  Figure  46  and  all  three  modes  appear  to 
be  the  second  bending  mode.  Since  two  fins  are  being  tested  as  one  structure, 
two  second  bending  modes  should  be  expected.  The  presence  of  three  modes  In 
the  frequency  band  of  294  to  328  H*  can  be  explained  by  slight  differences  In 
the  fins  which  were  of  different  vintage  and  probably  different  Internal 
structures.  Also,  small  differences  existed  In  the  short  (nominal  .191  Inch) 
cantilever  rods  for  each  fin.  The  two  frequencies  of  294  and  328  Hz  were 
selected  for  comparison  with  NASTRAN  predictions.  This  comparison  is  shown  in 
Figure  47.  Again,  acceptable  correlation  between  measured  and  predicted  modal 
frequencies  and  shapes  was  achieved. 
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The  last  Y-Z  direction  aode  to  be  aeasured  with  this  technique  and 
with  the  fins  in  this  orientation  was  located  at  482  Rs.  The  tabulated  data 
is  as  follows: 


TABLE  16.  Tin  Bending  Mode 

FREQUENCY  -  482  Rx 


LOG 

AMP 

REF  AMP 

PRASE 

C0S(Si-9?) 

SHAPE 

l(Clp) 

13.87 

2.70 

>  39.8 

+1.00 

+1.00 

2 

1.04 

1.60 

-  54.3 

O 

a 

1 

-  .01 

3 

13.35 

2.30 

+220.6 

-1.00 

-  .96 

4 

12.75 

2.07 

+227.7 

-  .99 

-  .91 

5 

2.57 

1.88 

-  18.5 

+  .52 

+  .10 

6 (base) 

1.30 

1.40 

+  23.3 

+  .95 

+  .09 

A  plot  of  this  aode  shape  is  shown  in  Figure  48.  The  aode  appesrs  to 
resemble  a  second  bending  aode  except  for  the  relativelY  low  vibration  levels 
near  the  base  of  the  fin.  Due  to  the  separation  of  frequency  froa  the  pre¬ 
viously  described  second  bending  aodes  in  the  300  Hz  frequency  range  and  the 
low  base  node  aaplitudes,  this  aode  was  classified  as  a  localized  effect  pro¬ 
bably  due  to  the  short  rod  section  used  to  cantilever  the  fins. 

B.  Fin  Random  Response  Y-Z  Direction. 

The  fins  were  retested  in  the  saae  direction  except  the  rsndoa 
response  technique  was  used  to  identify  the  aodes  of  vibri>tioo.  The  frequency 
response  is  shown  in  Figure  49.  The  frequencies  readily  identified  were  42.5, 
305,  326,  and  494  Rz.  Additional  acceleration  locations  were  used  during  these 
measurements.  As  shown  in  the  first  randoa  aode  plot.  Figure  50,  aeasureaents 
were  made  at  11  locations  with  position  1  being  the  fin  tip  and  position  11 
being  the  top  of  the  .650  inch  cube  block  used  to  cantilever  the  fins.  This 
plot  again  clearly  showed  the  first  aode  in  the  weak  axis  to  be  46  Hz.  The 
NASTRAN  prediction  was  not  repeated  since  the  correlation  was  excellent  when 
compared  to  the  previously  described  normal  aode  data. 

The  aode  shapes  for  the  300  Hz  frequency  range  resonances  are  shown 
in  Figures  51  and  52.  The  resonances  located  at  306  and  326  Hz  are  the  second 
bending  aodes.  These  modal  frequencies  were  found  to  be  294  and  328  Hz  when 
measured  by  the  normal  aode  aethod  which  is  again  excellent  correlation  bet¬ 
ween  the  two  techniques.  Note,  21  aeasureaents  were  made  using  the  randoa 
response  approach  so  both  fins  could  be  studied  simultaneously. 

The  last  frequency  studied  during  this  series  of  tests  was  located  at 
490  Hz  (see  Figure  53)  which  was  tht  ^«me  resonance  studied  at  482  Hz  i..  the 
normal  mode  technique.  Only  the  right  fin  was  studied.  Exactly  the  same 
minimal  response  of  the  nodes  near  the  base  is  apparent  as  previously  discussed, 
and  again  this  resonance  appeared  to  be  associated  with  vibration  of  the  outer 
portion  of  the  fin  with  the  cantilever  rod  remaining  essentially  stationary.  As 
dlscxissed  in  the  missile  airframe  body  bending  analysis,  the  structure  ateapted 
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to  approach  the  mode  shape  of  the  closest  ppedoalnant  aode  of  vibration  even 
when  a  subset  of  the  structure  or  subcoaponent  was  In  a  resonant  condition. 

C.  Fin  Noraal  Mode  X  Direction. 

The  fins  were  rotated  90  degrees  to  peralt  the  identification  of  the 
longitudinal  or  X  direction  modes  of  vibration.  The  sine  sweep  plot  for  the 
right  fin  is  shown  in  Figure  54.  The  vertical  axis  is  the  tip  acceleration 
divided  by  the  acceleration  at  the  base  of  the  fin  iaaediately  to  the  right  of 
the  cantilever  rod.  One  extreae,  alaost  classical,  resonance  Is  indicated  at 
at  81  Hz.  Subsequent  dwell  tests  pinpointed  the  frequency  to  be  82  Hz  as  given 
in  Table  17  and  Figure  55.  The  flexibility  was  provided  by  the  short  rod  used 
to  cantilever  the  fins  with  the  body  of  the  fins  essentially  acting  as  a  rigid 
oiass . 


TABLE  17.  Fin  First  Bending  Mode  X-Direction 
FREqCJENCT  •  82  Hz 


LOC 

AMP 

REP  AMP 

PRASE 

C0S(oi-o?) 

SHAPE 

l(tlp) 

12.90 

1.06 

-130.0 

+1.00 

+1.00 

2 

10.20 

.78 

-133.0 

+1.00 

+  .79 

3 

7.60 

.50 

-117.0 

+  .97 

+  .57 

4 

5.30 

.38 

-  88.0 

+  .74 

+  .30 

5 

3.10 

.38 

-  82.0 

+  .67 

+  .16 

6 (base) 

1.20 

.38 

-  63.0 

+  .39 

+  .04 

D.  Fin  Randoa  Response  X-Direction. 

The  frequency  response  function  is  shown  in  Figure  56  with  the  one 
resonance  Indicated  at  83  Hz.  Subsequent  aode  aeasureaents  were  conducted  at 
82  Hz  and  are  given  in  Figure  57.  Again,  excellent  correlation  was  provided 
between  the  noraal  aode  and  randoa  response  aeasureaent  techniques. 


V.  NASTRAN  FINITE  ELEMENT  MODEL  SENSITIVITY  ANALYSIS. 

Discussion  of  the  NASTRAN  analysis  has  been  liaited  in  the  preceding 
Chapters  and  only  the  final  results  have  been  presented.  Actually,  nuaerous 
iterations  were  made  of  the  basic  beaa  aodel  prior  to,  during,  and  subsequent 
to  the  testing.  In  one  case,  retesting  both  statically  and  dynaalcally  was 
accoaplished  to  better  understand  the  results  and  to  iaprove  correlation. 

This  sensitivity  analysis  is  laportant  because  it  demonstrates  soae  peculiari¬ 
ties  of  the  aodal  extraction  code  used  in  this  version  of  NASTRAN  and  soae 
considerations  that  should  be  given  such  analysis  and  test  prograas  in  the 
future. 

A  sunmtary  is  given  in  Table  18.  Seven  NASTRAN  models  were  studied  not 
counting  the  special  models  used  to  identify  the  body  bending  modes  discussed 
in  detail  in  Chapter  II.  Model  1  was  the  basic  model  and  was  used  for  the 
comparisons  between  predicted  and  measured  modes  except  for  the  fins.  This 
model  is  considered  the  baseline.  Note,  the  first  wing  bending  modes  had 
three  frequencies  (modes)  and  not  the  expected  four  associated  with  four 
wings.  Four  modes  are  given  for  the  fln(8)  first  bending  mode  and  all  higher 
frequency  bending  modes  for  both  the  wing(s)  and  fln(s). 

Model  2  is  a  repeat  of  the  first  model,  but  the  lower  wing  was  removed  as 
accoaplished  in  the  actual  test  sequence.  The  number  of  first  bending  modes 
reduced  to  two  and  the  higher  wing  modes  dropped  to  the  appropriate  number  of 
three.  A  slight  shift  in  frequencies  was  noted. 

Model  3  was  the  same  as  Model  2  except  that  new  wing  stiffness  measure¬ 
ments  became  available  from  other  FOG-M  structural  analyses.  With  these  new 
section  properties,  the  predicted  weak  wing  axis  bending  modes  dropped 
slightly,  but  correlation  was  neither  enhanced  nor  degraded  substantially. 

In  Model  4,  extremely  poor  correlation  was  achieved  between  predicted  and 
measured  first  fin  longitudinal  bending  mode  (378  Hz  predicted  and  81  Hz 
measured).  The  earlier  models  did  not  account  for  the  section  properties  of 
the  short  cantilever  beam  used  to  support  the  fins.  These  beam  elements  were 
added  and  the  predicted  modal  frequency  dropped  from  378  to  170  Hz  but  still 
did  not  correlate  satisfactorily.  This  significant  difference  was  corrected 
with  the  use  of  an  additional  static  test  and  will  be  discussed  with  Model  7. 

Models  5  and  6  were  continuations  of  the  effort  to  study  the  peculiarities 
of  this  particular  NASTRAN  modal  extraction  routine.  As  shown,  two  wings 
yield  one  first  wing  bending  mode,  and  one  wing  results  in  no  first  wing 
bending  mode.  This  observation  is  extremely  valuable  for  future  studies.  The 
first  mode  of  any  appendage  is  the  most  important;  consequently,  care  must  be 
taken  in  the  future  when  using  this  subroutine. 

Model  7  had  two  major  changes  from  the  previous  NASTRAN  models.  First, 
the  number  of  nodes  used  to  model  the  wings  and  fins  was  increased  signifi¬ 
cantly  in  an  attempt  to  Improve  the  correlation  of  the  second  and  tnlrd  weak 
axis  wing  bending  modes  and  the  second  weak  axis  fin  bending  modes  of 
vibration.  The  wing  nodes  were  Increased  from  4  to  20  and  the  fin  nodes 
increased  from  2  to  10.  Some  Improvement  was  obtained  but  the  frequency 
correlations  did  not  substantially  Improve.  The  mode  shapes  did  correlate 
extremely  well  for  all  these  modes. 
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TABLE  18.  NASTRAN  FEM  Sensitivity  Analysis. 
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The  other  change  In  the  model  was  that  static  tests  were  conducted  to 
measure  the  weak  and  strong  fin  axis  stiffness.  The  stiffness  was  translated 
Into  section  properties,  El  terms,  and  added  to  the  NASTRAN  model.  Marked 
Improvement  was  obtained  In  the  fin  longitudinal  (X~  direction)  frequency  pre¬ 
diction.  The  predicted  frequency  dropped  from  170  to  72  Hz  which  compares 
satisfactorily  with  the  81  Hz  measured.  The  new  weak  axis  stiffness  Increased 
the  first  bending  mode  prediction  from  41  Hz  (Model  4)  to  49  Hz  (Model  7)  com¬ 
pared  to  43  Hz  measured.  The  second  weak  axis  fin  bending  mode  prediction  was 
enhanced,  l.e.,  226  Hz  In  Model  4,  367  Hz  In  Model  7,  and  294-328  Hz  measured. 
Results  of  this  analysis  are  summarized  In  the  Lessons  Learned  Chapter. 
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VI.  LESSONS  LEARNED. 


A.  The  length  of  the  stinger  connecting  the  exciter  to  the  test  article 
should  be  kept  as  short  as  possible  (Fig.  43).  Longer  stingers  result  In 
relatively  low  frequency  lateral  vibrations  which  Induce  vibration  amplitudes 
Into  the  article  under  test. 

B.  The  procedures  for  using  the  force  transducer  recommended  by  the 
manufacturer  must  be  carefully  followed.  Extreme  dispersion  In  test  results 
was  found  when  the  tightening  torque  and  mating  lubricants  were  not  used  as 
specified.  A  factor  of  2  was  found  In  frequency  for  the  same  lower  frequency 
modes  of  vibration. 

C.  A  control  accelerometer  should  always  be  used  In  making  the  modal 
measurements  In  the  normal  mode  or  random  response  approach.  Limited  Instru¬ 
mentation  results  In  the  exciter  being  turned  on  and  off  during  the  process. 
Only  by  recording  a  reference  accelerometer  during  each  step  of  the  process 
can  a  baseline  be  established  to  give  validity  to  the  procedures. 

D.  A  major  problem  was  encountered  during  the  first  attempt  to  measure 
the  first  body  bending  mode  of  vibration.  The  data  was  Inconsistent  from  test 
setup  to  test  setup.  The  problem  was  Identified  as  loose  joints  which  hold 
the  cylindrical  body  shell  members  together.  There  are  seven  structural 
joints  which  connect  the  various  missile  subassemblies.  Screws  were  loose  and 
In  some  cases  missing.  A  concerted  effort  was  made  to  replace  defective  hard¬ 
ware  and  to  tighten  all  joints.  The  result  was  consistent  data  that  corre¬ 
lated  satisfactorily  with  the  NASTRAN  predictions. 

E.  Random  modal  testing  of  structures  Is  more  rapid  than  normal  modal 
testing  and  possibly  more  structured.  However,  a  definite  loss  for  the  feel  of 
the  data  can  result  since  a  tendency  Is  to  analyze  each  successive  resonance 
and  not  Identify  the  significance  of  the  modes  and  Interaction  of  structural 
parts  that  produce  less  significant  localized  vibrations.  A  definite  problem 
exists  In  the  random  method  when  a  subcomponent  Is  vibrating  at  a  frequency  In 
the  same  vicinity  as  the  structural  mode  under  Investigation.  Early  In  this 
test,  the  90  Hz  resonance  associated  with  the  cantilevered  fiber  optics  spool 
and  probably  with  the  glmbaled  TV  seeker  was  thought  to  be  the  first  bending 
mode  of  the  airframe.  Random  vibration  measurements  along  the  body  of  the 
missile  Indicated  the  90  Hz  to  be  the  first  bending  mode  of  the  airframe  which 
was  a  totally  erroneous  conclusion.  Typically,  when  a  subcomponent  Is 
vibrating  close  to  a  major  structural  mode,  the  mode  shape  will  assume  the 
shape  closest  to  the  nearest  major  structural  mode.  Again,  a  methodical 
approach  should  be  followed  by  starting  with  the  lowest  frequency  modes, 
correlating  these  modes  with  a  mathematical  model  one  at  a  time.  Investigating 
localized  resonances,  and  progressing  to  the  next  higher  frequency.  In  addi¬ 
tion  to  this  structured  approach.  It  Is  recommended  that  a  mixture  of  the  nor¬ 
mal  mode  and  random  response  techniques  be  used  to  cross  check  for  error  In 
future  tests. 
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F.  Great  care  should  be  exercised  in  using  some  of  NASTRAN's  modal 
extraction  subroutines.  As  discussed  In  the  previous  paragraph,  predictions 
for  the  critical  first  bending  modes  of  appendages  such  as  wings  and  fins  can 
be  lost  if  such  care  is  not  used. 

G.  All  tools  available  should  be  used  by  the  analyst/tester  to  correlate 
the  mathematical  models  such  as  static  load/def lection  tests  and  retests,  and 
iterations  of  the  mathematical  model;  and  repeat  nodal  testing  if  required. 

As  discussed  in  detail  in  Chapter  II,  MISSILE  AIRFRAME  ANALYSIS,  careful 
examination  and  separate  computer  runs  are  sometimes  necessary  to  identify  the 
less  flexible  body  bending  modes  when  masked  by  the  significantly  more 
flexible  vibration  of  wings  and  fins. 
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VII.  CONCLUSIONS 


The  first  eight  modes  of  vibration  of  a  FOG-M  mock-up  airframe  were 
measured.  One  additional  mode  of  vibration  was  indicated  during  a  previous 
vibration  survey  and  is  included  for  completeness.  These  nine  modes  of  vibra¬ 
tion  and  the  predicted  NASTRAN  structural  modes  are  summarized  in  Table  19. 

The  first  bending  modes  of  the  primary  missile  components  (wings,  fins,  and 
missile  body)  were  successfully  correlated  as  to  frequency  and  mode  shapes  and 
the  correlation  was  excellent  for  some  modes,  i.e.,  first  wing  bending. 

The  higher  bending  mode  shapes  correlated  well  all  the  way  through 
the  third  bending  mode  of  the  wing.  The  modal  frequencies  for  the  higher 
modes  of  vibration  did  not  correlate.  Possible  sources  for  this  lack  of 
agreement  were  offered  so  that  future  mathematical  models  can  be  improved. 
Further  iterations  of  the  mathematical  model  could  have  been  pursued  until 
more  successful  higher  mode  frequency  correlations  were  obtained,  but  this 
process  was  not  followed  for  the  following  reasons: 

-  These  tests  were  primarily  intended  to  demonstrate  the  versatility 
and  power  of  modal  analysis  and  testing.  This  goal  was  achieved. 

-  The  structure  under  study  was  not  a  flight  version  of  the  FOG-M 

-  The  first,  bending  modes  of  vibration  that  were  correlated  are  by 
far  the  most  important  and  represent  those  modes  of  vibration  where  the  maxi¬ 
mum  deflection  and  hence  the  strain  and  stress  energy  are  concentrated. 

The  stated  objectives  of  this  program  were  met  and  the  validity  of 
the  NASTRAN  model  was  established.  A  systematic  method  was  established  that, 
if  followed  in  the  future,  will  ensure  that  modal  tests  identify  critical 
modes  of  vibrations  and  establish  the  accuracy  of  NASTRAN  mathematical  models. 


TABLE  19.  Summary  FOG-M  Mocic-Up  Modal  Test/Analyals. 
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3rd  Wing  Bending  273  Not  506  Mode  Shape  Correlated 

Y-Z  Direction  330  Measured  Figure  23 


VIII.  RECOMMENDATIONS. 


A.  Repeat  this  process  for  an  lOE  FOG-M. 

B.  Model,  measure,  and  correlate  the  glmbaled  seeker  first  mode  of 
vibration  for  the  gimbaled  seeker. 

C.  Model,  measure,  and  correlate  the  fiber  spool  first  mode  of  vibra¬ 
tion  for  the  fiber  spool. 

D.  Apply  pneumatic  pressure  tc  the  fin  actuating  systems,  remove  any 
remaining  large  nonlinearity,  and  conduct  the  fin(s)  modal  tests  with  the  fins 
mounted  on  the  missile. 

E.  Measure  and  correlate  the  first  longitudinal  bending  mode  of  Che 
wingCs) . 

F.  Measure  Che  torsional  modes  of  vibration  of  the  wing(s)  and  fln(s). 
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Figure  3.  Transducer  locations  mlBslle  body  sine  swee 
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First  bending  nomal  mode  F  *  130 » 3  Hz 
first  bending  NASTRAN  F  -  147.1  Hz. 


O  NORMAL 
□  NflSTRRN 


iO.OQ  30.00  40.00  50.00 

INCH 


Second  bending  normal  mode  F  ■  267.0  Hz 
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NORMAL  MODE  f=175  HZ 
NASTRAN  MODE  f=90  HZ  TO  120  HZ 


Figure  20.  Left  wing  measured  and  predicted  second  bending  mode. 


- NORMAL  MODE  f=192  HZ 

- NASTRAN  MODE  f=90  HZ  TO  120  HZ 


Figure  21.  Left  wing  measured  and  predicted  second  bending  mode 


Figure  25.  Accelerometer  2,  random  vibration  15-210  Hz. 
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Figure  26.  Accelerometer  3>  random  vibration  15-210  Hz 
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Figure  28.  Acceleroaeter  5,  randoa  vlbratloa  13-210  Hz 
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RIGHT  WING 
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Figure  35.  Transducer  location  for  right  wing  test 
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Figure  43.  Fin  Y-Z  direction  cest  configuratioo. 
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Figure  44.  Fins  sine  sweep  frequency  range  10-510  Hz 
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Figure  49.  Fla  random  respoaae  fuactioo. 
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Figure  50.  Fin  random  response  first  bending  mode,  Frequency  "  46  Hz. 
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Mode  ♦  3  Damp;  .92  'A 

Freq:  305.81  Hz  View  ;<0,1,0> 


Figure  51.  Fin  random  response  second  bending  mode.  Frequency  ■  306  Hz. 


Figure  52.  Fin  random  response  second  bending  mode.  Frequency  ■  326  Hz. 
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Figure  54.  Fin  sine  sweep  X-direction. 


64 


MAG  TRHNS  *  1 


Figure  55.  F^ln  first  bending  mode  X-dlrectioo.  aormal  mode  method 
Frequeacy  ■  82  Hz.  ^  - - 
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Figure  56.  Fin  random  reapoose  function  X-direction. 
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